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Abstract 

Consider a polynomial F in m variables and a finite point ensemble 
S = Si X ■ ■ ■ X Sm . When given the leading monomial of F with respect to 
a lexicographic ordering we derive improved information on the possible 
number of zeros of F of multiplicity at least r from S. We then use this 
information to design a list decoding algorithm for a large class of affine 
variety codes. 

Keywords. Affine variety code, list decoding, multiplicity, Schwartz- 
Zippel bound. 

1 Introduction 

In this paper we study affine variety codes over Fq in the special case where the 
variety is S*! x • ■ • x Sm- Here, Si C Fq, i = 1, . . . ,m. More formally write 

S^SiX---xSm = {Pu...,PlS\} 

and consider the evaluation map 

evs :F,[Xi,...,X„,] ^fI^I, evs(F) = (F(Pi), . . . , F(F|5|)). 

Let 

MC{Xl^...X^- \0<z, < |5,|,j = l,...,m} 
and define the affine variety code 

£;(M, S) = SpanFjevs(M) | 7\f e M}. 

Clearly, the generalized Reed-MuUer codes is a special example of such a code. 
Other well-known examples are the hyperbolic codes [5], Reed-Solomon prod- 
uct codes [T] , (generalized) toric codes [TUl [HI [3] , and Joyner codes [Tl] . In the 
present work we show that the dimension of E{M., S) equals |M| and we present 
a lower bound on the minimum distance which turns out to be sharp when |M| 
satisfies certain reasonable criteria. We then present a decoding algorithm for 
E(Wl, S). The algorithm is a straightforward generalization of the Guruswami- 
Sudan decoding algorithm for Reed-Solomon codes |i9l, except that it involves 
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a preparation step. A main ingredient in the preparation step is information 
about how many zeros of multiplicity at least r a polynomial F{Xi, . . . ,Xm) 
can have given information about its leading monomial with respect to a lexi- 
cographic ordering. 

It is well-known that the generalized Reed-Muller codes can be viewed as sub- 
field subcodes of Reed-Solomon codes. As demonstrated by Pellikaan and Wu 
[m [13] this observation leads to an efficient decoding algorithm of general- 
ized Reed-Muller codes via the Guruswami- Sudan decoding algorithm for Reed- 
Solomon codes. When the variety does not equal S" = x ■ • • x F^, or when 
M is not chosen as all the monomials of degree up to some number, such a neat 
result does not hold any more. 

Pellikaan and Wu also presented a direct interpretation of the Guruswami-Sudan 
decoding algorithm for generalized Reed-Muller codes. Even though Augot et 
al. improved the analysis of the latter algorithm dramatically [U [5] the algo- 
rithm that uses the subfield subcode is still superior when decoding generalized 
Reed-Muller codes. The analysis by Augot et al. uses a generalization of the 
Schwartz-Zippel bound to also deal with multiplicity. A proof of this bound 
was given by Dvir et al. in [5]. Augot et al. then used their insight to general- 
ize the second algorithm by Pellikaan and Wu (the direct interpretation of the 
Guruwami-Sudan algorithm) to Reed-Muller like codes over 5' = S'l x • ■ • x Sm, 
when Si — ■ ■ ■ — S*™, and to Reed-Solomon product codes as well. Based on 
the generalized Schwartz-Zippel bound they estimated the decoding radius. 
In the present work we improve upon the methods from [5] to derive more de- 
tailed information on how many zeros of prescribed multiplicity a polynomial 
can have. Rather than using only information on the degree of the polynomial 
we use information on the leading monomial with respect to a lexicographic 
ordering. We present both a recursive algorithm to find such bounds and also 
closed formula expressions for the case of polynomials in two variables. 
The interpolation polynomial Q{Xi, . . . , Xm, Z) in the list decoding algorithm 
of the present paper can be viewed as a polynomial in Z with coefhcients from 
Fg[Xi, . . . ,Xm]- Having fixed the code to be used, in a preparation step we 
determine sets from which the coefficients must be taken. Here, we use our 
improved insight on how many zeros with prescribed multiplicity a polynomial 
can have given information about its leading monomial. The idea of a prepa- 
ration step comes from [7J where an interpretation of the Guruswami-Sudan 
decoding algorithm without multiplicity was described for order domain codes. 
Our experiments show that the method of the present work is an improvement 
upon the situation where only the generalized Schwartz-Zippel bound is used 
for the design and the analysis of the decoding algorithm. Our algorithm works 
for codes of not too high dimensions. For small dimensions the algorithm often 
decodes more than d/2 errors. 

2 Parameters of the codes 

Throughout the paper we shall use the notation Si = \Si \ for i = 1, . . . , m. Also 
we shall write n = \S\ as clearly the length of i?(M, S) equals [S'l. First we show 
how to find the dimension of the code. 

Proposition 1. The dimension of E(M., S) equals |M|. 
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Proof. We only need to show that 

{evs(Xi\...,X;,"0 |0<i, <s„j = l,...,m} 

constitutes a basis for FJ^' as a vectorspace over Fq. For this purpose it is 
sufficient to show that the restriction of ev to 

{G(Xi, . . . , X,„) I deg^^ (G) < s„ z = 1, . . . , m} (f ) 

is surjective. Given (ai, . . . , a„) £ F^ let 

- a 



Here, we have used the notation P„ = (P^J-^'', . . . , Pi^-*), v = 1, . . . ,n. It is 
clear that evs(P) = (ai, . . . ,a„) and therefore evs ■ Fq[Xi, . . . ,X,n] F^ is 
surjective. Consider an arbitrary monomial ordering. Let R{Xi, . . . , Xm) be 
the remainder of F{Xi, . . . , X^) after division with 

. . . ,X„i), ■ ■ . ,A,n{Xi, . . . ,Xm)}- 

Here, A,{Xi, . . . = UaeS^^^-"'): * = 1, • ■ • , Clearly, P(P,) = R{P,) = 

ai, i = I, . . . ,n. Hence, the restriction of evs to ^ is indeed surjective. □ 

We next show how to estimate the minimum distance of E{M.,S). The 
Schwartz-Zippel bound [151 113 E] is as follows: 

Theorem 2. Given a lexicographic ordering let the leading monomial oj F{Xi, . . . , Xm) 
he X^^ ■ ■ ■ • The number of elements in S ^ Si x ■ ■ ■ x Sm that are zeros of 
F is at most equal to 

ilS2 ■■■Sm+ Sl?2S3 ' ' • Sm H h «! • • ■ Sm^iim- 

The proof of this result is purely combinatorial. Using the inclusion-exclusion 
principle it can actually be strengthened to the following result which is a special 
case of the footprint bound from Grobner basis theory: 

Theorem 3. Given a lexicographic ordering let the leading monomial of F{Xi, . . . , Xm) 
be Xl^ ■ ■ ■ Xl^ . The number of elements in S = Si x ■ ■ ■ x Sm that are zeros of 
F is at most equal to 

n - (Si - n)(s2 - ^2) • • • [Sm - im)- 

Proposition 4. The minimum distance of E{Wl, S) is at least 

min{(si - h)(s2 - ^2) ■ • • {sm - irn)\Xl' ■■■X'^ M}. 

If for every M e M it holds that N e M for all N such that N\M then the 
bound is sharp. 
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Proof. The first part follows from Theorem [31 To see the last part write for 
i = I, . . . ,m, Si = {b'"i \ . . . , &|5. [}• The polynomial 

F(xi,...x„)=nn(^"-^r) 

11=1 j=i 

has leading monomial Xl^ ■ ■ ■ with respect to any monomial ordering and 
evaluates to zero in exactly n — (si — ii)(s2 — 12) ■ ■ ■ (sm — im) points from S. 
Finally, any monomial that occurs in the support of F is a factor of Xl^ ■ ■ ■ X^ . 

□ 



3 Bounding the number of zeros of multiplicity r 

In the following let X = {Xi, X^) and f = (Ti, . . . , T„). 

Definition 5. Let F be any field. Given F{X) e F[X] and k e NJ^ then 
the k'th Hasse derivative of F, denoted by F'^''^{X), is the coefficient of in 
F{X + T). In other words 

F[X + f) = ^F('=)(X)f^. 

The concept of multiplicity for univariate polynomials is generalized to mul- 
tivariate polynomials in the following way: 

Definition 6. For F{X) e F[X]\{0} and a £ F'" we define the multiplicity 
of F at a denoted by mult{F, a) as follows: Let r be an integer such that for 
every k — {ki, . . . ,km) G N™ with ki + ■ ■ ■ + km < r, F^'^^S) = holds, but for 
some k = (fci, . . . , km) £ N™ with ki + ■ ■ ■ + km = r, i^('^^(a) ^ holds, then 
mult{F, a) = r. If F — then we define mult{F, a) — 00. 

Elaborating on the results in [5 we find: 

Theorem 7. Let F{X) G F[X] be a non-zero polynomial and let X^^ ■ ■ ■ X^j^ 
be its leading monomial with respect to a lexicographic ordering. Then for any 
finite sets Si, ... , Sm ^ F 

mult{F, a) < iiS2 • • • Sm + Sl«2S3 • • • H h si • • ■ Sm-lim- 

aGSi X ••• X Sm 

Theorem[2]now generalizes to the following result which we call the Schwartz- 
Zippel bound. 

Corollary 8. Let F{X) e F[X] be a non-zero polynomial and let X]^ ■ ■ ■ be 
its leading monomial with respect to a lexicographic ordering. A ssume Si, ... , Sm ^ 
F are finite sets. Then over Si x ■ ■ ■ x Sm the number of zeros of multiplicity 
at least r is less than or equal to 

(«lS2 • • • Sm + 51*253 • ' • Sm H h Si • • • Sm-lim) /r. (2) 
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Just as we were able to improve upon Theorem [2] by using the inclusion- 
exclusion principle we will also be able to improve upon Corollary [S] However, 
now the situation is much more complex and therefore the inclusion-exclusion 
principle is no longer sufficient. What we need to strengthen Corollary [5] is the 
following rather technical function: 

Definition 9. Let r e N, ii, . . . , G Nq. Define 

D{ii,r,si) = min{[^J,si} 

and for m > 2 

max <^ (s,„ - Ul Ur)D{ii, . . r,si, . . .,Sm-i) 

+ uiD{ii,... ,im-i,r - 1, si, . . . , s,„-i) H 

-I- Ur-lD{ii, . . . ,im-l, l,Sl, . . . , Sm-l) "I- U^Sl • • • S,„_i 

where 

{{ui, . . . , Ur) e Nq \ Ul + ■ ■ ■ + Ur < s„i and Ui + 2U2 + • • • + rUr < im}- 

Theorem 10. For a polynomial F{X) G F[X] let ■ ■ ■ X^^'^ be its leading 
monomial with respect to the lexicographic ordering with Xm -<■■■< Xi . Then 
F has at most D{ii, . . . , im, r, si, . . . , Sm) zeros of multiplicity at least r in Si x 
• • • X Sm ■ The corresponding recursive algorithm produces a number that is at 
most equal to the number found in Corollary\^ and at most equal to Si ■ ■ ■ Sm- 

Proof. The proof involves a modification of the method described in Due 
to lack of space we do not include it here. □ 

Remark 11. Given (ii, . . . , im, r,si, . . . , Sm) with [ii/sij + • ■ • + [im/sml > r 
then there exist polynomials with the leading monomial being Xl^ ■ ■ ■ such 
that all points in Si x ■ ■ ■ x Sm o.fe zeros of multiplicity at least r. Hence, we 
need only apply the algorithm to tuples (ii, . . . ,im) such that 

Ui/siJ H 1- [im/Sm\ < r. (3) 

In a number of experiments listed in we calculated the value D(ii, . . . , i„i, r, . . . , g) 
for various choices of m, q and r and for all values of (ii, . . . , im) satisfying ([3]). 
We here list the maximal attained improvement obtained by using Proposi- 
tion [TU] rather than using Corollary [51 We do this relatively to the number of 
points q™. In other words we list in Tabled] the value 

( max {min{(ii + ■■■ im)q"'~^/r, (?™} - D{ii, . . . ,im,r,q, . . . ,q)}] /q™ 
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Table 1; Maximum improvements relative to q™; truncated 



m 2 3 4 
r2345234523 

2 0.25 0.25 0.25 0.25 0.25 0.375 0.375 0.375 0.312 0.375 

3 0.222 0.222 0.222 0.222 0.296 0.296 0.296 0.296 0.296 0.333 

4 0.187 0.187 0.187 0.187 0.281 0.25 0.25 0.265 0.316 0.289 
q 5 0.24 0.16 0.16 0.2 0.256 0.256 0.232 0.24 0.307 0.288 

7 0.204 0.204 0.163 0.142 0.279 0.244 0.227 0.209 0.299 0.276 

8 0.234 0.203 0.171 0.140 0.275 0.25 0.214 0.203 0.299 0.275 



for various choices oi m,q,r. 

The experiments also show distinct average improvement. This is illustrated in 
Table [5] where for fixed q,r,m we list the mean value of 

min{(zi + h im)q"'^^ , q"'} - Djii, . . . ,im,r, q, . . . , q) 

min{(ii + --- + i„Oq"-i,g'"} 

The average is taken over the set of exponents (ii, . . . , im) where [ii/gj + 
1- [im/ql < r holds. 



Table 2: The mean value of ([4]); truncated 
m 2 3 4 

r2345234523 

2 0.363 0.273 0.337 0.291 0.301 0.300 0.342 0.307 0.248 0.260 

3 0.217 0.286 0.228 0.236 0.194 0.224 0.213 0.214 0.158 0.177 

4 0.191 0.197 0.232 0.195 0.158 0.169 0.180 0.172 0.125 0.135 
q 5 0.155 0.167 0.174 0.197 0.139 0.145 0.148 0.153 0.110 0.116 

7 0.128 0.137 0.138 0.138 0.119 0.122 0.121 0.119 0.093 0.098 

8 0.126 0.127 0.134 0.126 0.114 0.115 0.113 0.111 0.089 0.093 



The values D{ii, . . . , im, si, . . . , Sm) may sometimes be time consuming to 
calculate. Therefore it is relevant to have some closed formula estimates of these 
numbers. We next present such estimates for the case of two variables. By Re- 
mark [11] the following proposition covers all non-trivial cases: 

Proposition 12. For fc = 1, . . . , r — 1, D{ii,i2, r, si, S2) is upper bounded by 

(C.l) S2— + r 

r r r — k 

if [r — A;)^q-j-si < ii < [r — k)si and < i2 < ks2, 

(C.2) S2- + ((fc + 1).S2 - 12){^ --) + {l2 - kS2){si ~ ^) 

r r — k r r 

if (r — A;)^-j-si < *i < (t" — fc)si and ks2 < 12 < (fc + l)s2j 

(C.3) s2^ + ^is~-'-^) 
r k + 1 r 

if [r — k — l)si < ii < {r — fc)^q-j-si and < i2 < {k + l)s2- 

Finally, 
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(C.4) D{ii,i2,r,si,S2) = S2L— J +^2(31 - [—J) 

r r 

if Si{r — 1) < ii < Sir and < 12 < S2- 
The above numbers are at most equal to min{(iiS2 + Sii2)/r, 8182}- 

Proof. The estimates are found by treating ii and 12 as rational numbers rather 
than integers. Due to lack of space we do not give the details here. □ 

4 The decoding algorithm 

The main ingredient of the decoding algorithm is to find an interpolation poly- 
nomial 

Q{Xi, . . . , X,n, Z) = Qo{Xi, . . . , Xm)+Qi{Xi, . . . , X„i)Z+- ■ ■+Qt{Xi, . . . , X,„)Z* 

such that Q{Xi, . . . , Xm, F{Xi, . . . , Xm)) cannot have more than n — E different 
zeros of multiplicity at least r whenever Supp(F) C M. The integer E above is 
the number of errors to be corrected by our list decoding algorithm. To fulfill this 
requirement we will define appropriate sets of monomials B{i, E,r), i — 1, . . . ,t 
and then require Qi{Xi, . . . , Xm) to be chosen such that Supp((5i) C B{i, E, r). 
Rather than using the results from the previous section on all possible choices 
of F{Xi, . . . ,Xm) with Supp(F) C M we only consider the worst cases where 
the leading monomial of F is contained in the following set: 
Definition 13. 

M = {Af e M I if N and M\N then M = N}. 

Hence, M is so to speak the boarder of M. 

Definition 14. Given positive integers i,E,r with E < n let 

B{i, E, r) = {/V e A(r, m) | Dr{KA'P) <n- E for all M e M}. 

Here, Dr{Xl^ , • • • , X^^) can either be D{ii, . . . , i™, r, si, . . . , Sm) or in the case 
of two variables it can be the numbers from Proposition \1'A Another option 
would be to let Dr{Xl^ , . . . , X^^) be the number in With a reference to 

Remark \ll\ we have defined 

A(r, m) = {X\' • • • I Ui/siJ + • • • + [im/smj < r}. 
The decoding algorithm calls for positive integers t, E, r such that 
t 

J2\B{i,E,r)\ > nN{m,r), (5) 
1=1 

where N{m, r) = is the number of linear equations to be satisfied for a 

point in F™+^ to be a zero of Q{Xi, . . . , Xm, Z) of multiplicity at least r. As 
we will see condition ^ ensures that we can correct E errors. We will say that 
{t,E,r) satisfies the initial condition if given the pair {E,r), t is the smallest 
integer such that ([SJ is satisfied. Whenever this is the case we define B'{t, E, r) 
to be any subset of B{t, E, r) such that 

t-i 

\B{i,E,r)\ + \B'{t,E,r)\ = nN{m,r) + 1. 

i=l 

Replacing B{t, E, r) with B'{t, E, r) will lower the run time of the algorithm. 
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Algorithm 15. Input: 

Received word r = (ri, . . . , r„) £ F^. 

Set of integers {t, E, r) that satisfies the initial condition. 

Corresponding sets B{1, E,r) . . . , B{t — 1, E, r), B'(t, E, r). 

Step 1 

Find non-zero polynomial 

Q{Xi, . . . , X„iZ) = Qo{Xi, . . . , Xm)+Ql{Xi, . . . , Xm)Z+- ■ ■+Qt{Xi, . . . , Xm)Z* 

such that 

1. Supp{Q^) C B{i,E,r) fori = l,...,t-l and Supp{Qt) C B'{t,E,r), 

2. {Pi,ri) is a zero of Q{Xi^ . . . , X„n Z) of multiplicity at least r for i — 
1, . . . 

Step 2 

Find all F{Xi, . . . , Xm) G Fq[Xi, . . . , X,n] such that 

[Z - . . . , Xra))\QiXi, ...,Xra, Z). (6) 

Output: 

A list containing (F{Pi), . . . , F{Pn)) for all F satisfying (0). 

Theorem 16. The output of A laorithm \15\ contains all words in E(Nl, S) within 
distance E from the received word r. Once the preparation step has been per- 
formed the algorithm runs in time 0{n^) where n — "-("^x)' ff*'^^'^ mul- 
tiplicity r the maximal number of correctable errors E and the corresponding sets 
B{1, E,r),..., B{t-l, E, r), B'{t, E, r) can be found in time 0(nlog(n)r™s'|M|/cr) 
assuming that the values of the function Dr are known. Here a = max{ii + 
• ■ • + im\Xl^ ■ ■ ■ X^ G M} and s' — niax{si, . . . , Sm}- 

Proof. The interpolation problem corresponds to n homogeneous linear equa- 
tions in n + 1 unknowns. Hence, indeed a suitable Q can be found in time 
Now assume Supp(F) C M and that distff(evs(J^), f) < E. Then 
Pj is a zero of Q{Xi, . . . , Xm, F{Xi, . . . , Xm)) of multiplicity at least r for 
at least n ~ E choices of j. By the definition of B{i,E,r) this can, how- 
ever, only be the case if Q{Xi, . . . , Xm, F{Xi, . . . , X„i)) ~ 0. Therefore, Z — 
F{Xi, . . . , Xjn) is a factor in Q{Xi, . . . , Xm, Z). Finding linear factors of poly- 
nomials in {Fq[Xi, . . . ,Xm])[Z] can be done in time 0{n^) by applying Wu's 
algorithm in [IS] (see [H P- 21]). □ 

5 Examples 

Example 17. In this example we consider a point ensemble S — Si x S2 with 
si = 128 and S2 = 64. We consider codes E{m,S) where M = {X^^X*" | 
ii -\- 2^2 < u}. Note, the weight 2 which corresponds to the number Si/s2- The 
performance of A laorithm \15\is independent of the field in which Si and S2 live. 
In Table we list the number of errors that we can correct when the function 
^r(ii,'*2) - is chosen to 6e 12, r, 80, 80) (column D), - is the closed formula 
expression from Proposition]!^ (column C), - or is the Schwartz- Zippel bound 
(column S), respectively. The row \ {d— 1)/2J corresponds to half the minimum 
distance and the row Dim. is the dimension of the code. 
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Table 3: Error correction capabilities for the codes in Example [T71 



u 


3 


4 


7 


20 


r 


DCS 


DCS 


DCS 


DCS 


2 


5129 5105 4895 


4799 4777 4575 


4143 4124 3871 


2487 2475 2175 


3 


5367 5333 5205 


5048 5016 4906 


4407 4381 4245 


2855 2833 2666 


4 


5474 5438 5343 


5180 5143 5071 


4566 4535 4431 


3060 3031 2927 


9 


5653 5617 


5390 5361 


4817 4785 


3415 3384 


20 


5757 5740 


5509 5494 


4959 4943 


3609 3599 


L^J 


3999 


3967 


3871 


3455 


Dim. 


6 


9 


20 


121 


Table 4: Error correction capabilities 


for the codes in 


Example [T51 


iki,k2) 


(4,7) 


(5,9) 


(8,15) 


(21,41) 


r 


DCS 


DCS 


DCS 


DCS 


2 


4036 4015 3519 


3655 3639 3071 


2820 2808 2111 


1061 1055 


3 


4289 4261 3903 


3911 3885 3498 


3077 3058 2602 


1183 1171 533 


4 


4411 4381 4111 


4042 4011 3727 


3214 3187 2847 


1310 1291 831 


9 


4598 4487 


4244 4124 


3455 3313 


1567 1365 


20 


4711 4662 


4364 4310 


3588 3526 


1704 1615 




3720 


3599 


3248 


1935 


Dim. 


28 


45 


120 


861 



Example 18. This is a continuation of Example [7^ We consider the same 
point ensemble but choose M = {Xl^X2^ | ii < fci,i2 < /c2}. That is, we 
consider Reed-Solomon product codes. In Table j^] we list the number of errors 
that can be corrected by Alaorithm \15\ 

Example 19. In this example we consider a point ensemble S — Si x S2 with 
si = S2 = 80. We consider codes E{M,S) where M = {X{^Xl^ Mi + 12 < 
w} for various values of u. That is, we consider Reed-Muller like codes. The 
performance of A laorithm \15\ is independent on the field in which Si and S2 live. 
In Table [5| we list the number of errors that we can correct when the function 
Dr{ii,i2) ^ is chosen to &e Z?(ii, 12, r, 80, 80) (column D), - is the closed formula 
expression from Provosition (column C), - or is the Schwartz- Zippel bound 
(column S), respectively. Column A corresponds to a bound from on what 
can be achieved by applying their algorithm. The row A^o corresponds to what 
could theoretically be achieved by the algorithm in '21 if one uses high enough 
multiplicity. 

Assuming Si, S2 ^ F128, i?(M, S) corresponds to a puncturing of the generalized 
Reed-Muller code RMi2s{u,2). This suggest that as an alternative to using 
Alaorithm \15\ one could decode with respect to RMi2s{u, 2) treating the punctured 
points as errors. The best known algorithm to decode RMi2%{u, 2) for the values 
of u considered in this example is the algorithm by Pellikaan and Wu which 
uses the sub field subcode approach. Row PW of the table explains what can be 
achieved by this alternative approach. 
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Table 5: Error correction capabilities for the codes in Example [T^ 



u 


3 


4 


7 




20 




r 


D C S A 


D C S A 


D G S A 


D 


C S 


A 


2 


3594 3571 3399 3310 


3317 3297 3119 2999 


2693 2679 2479 2302 


1279 


1279 999 


585 


3 


3791 3765 3679 3604 


3524 3499 3413 3323 


2943 2918 2799 2692 


1575 


1559 1439 


1138 


4 


3899 3869 3799 3758 


3647 3618 3559 3492 


3080 3058 2979 2896 




1728 1639 


1428 


9 


4072 4053 4027 


3837 3813 3792 


3315 3297 3253 




2053 2035 


1935 


20 


4171 4163 4152 


3946 3939 3926 


3444 3435 3418 




2219 2211 


2169 




4257 


4042 


3558 




2368 




PW 


3891 


3503 


2568 











3079 


3039 


2919 




2399 




Dim. 


10 


15 


36 




231 
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